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ABSTRACT 


Measurements  are  presented  of  boundary  layers  with  embedded  vortices  and  with 
film  cooling  for  freestream  velocities  of  15,  and  11  m  s.  Measurements  of  a  boundary- 
layer  with  embedded  vortex  and  and  no  film  cooling,  and  of  a  boundary  layer  with  film 
cooling  but  no  vortex  are  presented  for  freestream  velocity  of  15  m  s.  Plots  of  total 
velocity.  V,  streamwise  velocity,  Vv  secondary  flow  vectors,  total  pressure.  PQ,  and 
streamwise  vorticity  are  presented  for  many  of  these  test  conditions. 

The  results  show  that  the  embedded  vortices  completely  dominate  the  flow  field 
in  boundary-  layers  with  film  cooling.  This  is  indicated  from  the  plots  of  V,  V  and  PQ 
which  show  the  effects  of  film  cooling  to  be  completely  decimated  in  the  vicinity  of  the 
vortex. 

In  order  to  conduct  this  study,  a  five-hole  pressure  probe  was  calibrated  for  pitch 
and  yaw.  The  probe  was  then  used  to  measure  five  pressures  associated  with  the  flow. 
From  these  pressures,  total  velocity  and  the  x.  y.  and  z  components  of  velocity  were 
determined. 

A  boundary  layer  profile  was  conducted  to  verify  the  calibration  of  the  pressure 
probe,  measurement  procedures,  and  velocity  computations.  The  results  show  expected 
boundary  layer  behavior  with  a  small  \\.  and  V'z  component. 
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I.  INTRODUCTION 


The  increasing  need  for  greater  efficiency  in  gas  turbine  engines  has  resulted  in 
higher  turbine  inlet  temperatures.  Consequently,  combustor  liners  and  turbine  blading 
are  subjected  to  greater  amounts  of  thermal  stress,  thermal  fatigue,  and  creep.  At 
present,  gas  turbines,  such  as  those  associated  with  military  applications,  have  inlet 
temperatures  as  high  as  1800  -  2000  degrees  C  (3270  -  3632  F)  with  pressures  of  40 
atmospheres. 

Turbine  parts  may  be  protected  from  heat  loads  resulting  from  exposure  to  gas  at 
high  temperatures  by  using  a  coolant  within  turbine  passages  and  along  turbine 
surfaces.  Convection  cooling,  impingement  cooling,  transpiration  cooling,  and  film 
cooling  are  used  for  this  purpose.  Although  it  is  possible  for  the  cooling  medium  to  be 
a  substance  such  as  liquid  water  or  Freon  -  12,  most  gas  turbine  arrangements  utilize 
engine  air  bled  off  from  the  compressor  and  rerouted  to  the  turbine  nozzles  and 
blading.  With  film  cooling,  compressor  air  is  ejected  from  surfaces  of  blades  and 
vanes.  The  film  coolant  then  protects  metal  surfaces  by  forming  a  protective  insulating 
film  between  the  blades  and  the  mainstream,  and  by  acting  as  a  heat  sink. 

The  (low  through  a  turbine  cascade  is  extremely  complex.  Efforts  to  analytically 
model  the  flow  are  successful  only  when  the  fluid  is  considered  to  be  inviscid.  When 
the  effects  of  viscosity  are  included,  the  analysis  is  much  less  effective  [Ref.  1).  Because 
of  the  difficulty  involved  in  analytical  representation  of  the  viscous  portions  of  cascade 
flows,  much  of  the  work  in  this  area  is  experimental  in  nature. 

Flow  visualization  studies  by  Fferzig,  et  al.,  [Ref.  2]  were  among  the  first  to  show 
the  complexity  of  flow  within  the  cascade.  Flow  visualization  studies  by  Langston,  et 
al..  [Ref.  3]  and  Marchal  and  Sieverding,  [Ref.  4[  also  show  the  detailed  development 
and  complexity  of  the  flow  through  the  turbine  cascade.  More  recently,  Sieverding  and 
Van  Den  Bosche,  [Ref.  5]  have  used  color  smoke-visualization  to  study  the  evolution  of 
(low  in  cascades. 

Figure  1.1  from  [Ref.  6]  shows  the  various  secondary  flows  associated  with  flow 
in  turbine  cascades.  As  the  inlet  boundary  layer  approaches  the  blade  s  leading  edge,  a 
horseshoe  vortex  is  formed.  The  point  at  which  this  formation  occurs  is  the  saddle 
point  which  is  clearly  shown  in  Figure  1.2  from  [Ref.  4],  One  leg  of  the  horseshoe 


vortex  moves  initially  near  the  side  of  the  blade,  and  then  through  the  blade  passage 
where  it  is  referred  to  as  the  passage  vortex.  In  cascade  flow,  the  passage  vortex  is 
composed  of  1)  fluid  from  the  pressure  side  leg  of  the  horseshoe  vortex,  2)  the 
crossflow  from  the  endwall  boundary  layer,  and  3)  entrained  fluid  from  the  mainstream 
flow,  [Ref.  7].  As  the  passage  vortex  continues  through  the  cascade,  it  is  forced  by  the 
pressure  gradient  to  the  opposite  side  of  the  passage  near  the  suction  side  of  the 
adjacent  blade.  This  shift  is  clearly  shown  in  Figure  1.3  from  [Ref.  2j.  The  second  leg 
of  the  horseshoe  vortex  follows  the  suction  side  of  the  blade  and  has  a  sense  of 
rotation  opposite  to  that  of  the  pressure  side  vortex.  This  vortex  moves  away  from  the 
corner  and  is  generally  believed  to  be  smaller  in  size  than  the  passage  vortex.  As  can 
be  seen  from  Figure  1.4  from  [Ref.  8]  the  passage  vortex  makes  approximately  one 
rotation  as  it  passes  through  the  cascade.  In  Figure  1.1  the  number  of  revolutions  has 
been  exaggerated  for  clarity. 

The  objective  of  this  thesis  is  to  study  the  effects  of  embedded  vortices  on  a  film 
cooled  turbulent  boundary  layer.  In  order  to  understand  the  effects  that  this  complex 
flow  field  has  upon  heat  transfer,  vortex  characteristics  and  their  interaction  with 
surrounding  flow  must  be  understood.  In  this  study,  a  five-hole  pressure  probe  is  used 
to  measure  vortex  characteristics. 

Extensive  procedures  for  calibration  and  qualification  of  the  five-hole  pressure 
probe,  and  its  use  in  measuring  three-dimensional  flows  are  first  discussed.  Results  of 
six  different  rest  are  then  given:  1)  a  baseline  measurement  in  a  developing  boundary 
layer,  2)  a  boundary  layer  with  film  cooling  only  at  a  freestream  velocity  of  15m  s.  3)  a 
boundary  layer  with  a  single  embedded  vortex  and  without  film  cooling  at  a  freestream 
velocity  of  15  m,'s,  4)  a  boundary  layer  with  embedded  vortex  and  film  cooling  at  a 
freestream  velocity  of  15  m/s,  5)  a  boundary  layer  with  an  embedded  vortex  and  no 
film  cooling  at  a  freestream  velocity  of  20  m/s,  and  6)  a  boundary  layer  with  embedded 
vortex  and  film  cooling  at  a  freestream  velocity  of  11  m  s.  Cases  3.  and  4  provide  the 
most  informative  results  of  the  study;  cases  5  and  6  are  incidental. 


II.  EXPERIMENTAL  APPARATUS 


A.  WIND  TUNNEL 

The  wind  tunnel  pictured  in  Figure  2.1  is  now  described.  It  is  an  open-circuit 
blower  tunnel  used  to  provide  uniform  flow  at  the  nozzle  exit. 

1.  Description 

The  wind  tunnel  is  designated  the  NPS  Shear  Layer  Research  Facility  (SLRF) 
and  was  built  by  Aerolab.  It  is  designed  to  provide  uniform  flow  with  a  minimum 
amount  of  turbulence  intensity.  It  is  designed  with  numerous  pressure  taps  and  four 
38  x  20.3  cm  (15  x  8  in.)  access  ports  along  each  of  the  of  the  side  walls.  The  height  of 
the  top  wall  is  adjustable  to  permit  changes  in  the  pressure  gradient  along  the  length  of 
the  test  section.  Additionally,  the  top  wall  contains  numerous  instrument  ports  for  the 
measurement  of  various  flow  characteristics. 

The  air  speed  through  the  tunnel  can  be  adjusted  from  5  m/s  to  40  m  s.  The 
blower  exit  slips  into  the  inlet  end  of  the  wide-angle  diffuser  with  1.6  mm  of  clearance 
so  that  the  fan  is  isolated  from  the  body  of  the  wind  tunnel.  The  diffuser  inlet  contains 
a  filter  and  nozzle.  The  test  section  is  3.048  m  (10  ft.)  long  and  0.6096  m  (2  ft.)  wide. 
The  top  is  fabricated  from  Lexan  sheet  (4.76  nun  thick),  continuously  sealed  with 
neoprene  along  the  edges.  The  tunnel's  bottom  wall  consists  of  one  1.2192  m  (4  ft.) 
long  and  three  0.6096  m  (2  ft.)  long  removeable  and  replaceable  sections.  These 
sections  are  all  0.6096  m  (2  ft.)  wide  and  are  sealed  with  "O"  rings  around  the  sides. 
Further  discussions  of  the  wind  tunnel  are  contained  in  [Ref.  9]  and  [Ref.  10:  p.  38]. 

2.  Qualification  and  Performance 

Extensive  qualification  test  of  the  Shear  Layer  Research  Facility  were 
conducted  by  Ligrani,  [Ref.  11].  Results  show  that  the  variation  of  total  pressure  at 
the  exit  plane  of  the  nozzle  is  less  than  0.4%  at  26  m  s  and  34  m/s.  Mean  velocity 
varies  less  than  0.7%  for  the  same  mean  freestream  speeds.  From  five-hole  pressure 
probe  measurements,  the  velocity  angle  deviation  is  nowhere  greater  than  about  0.6 
degrees  at  the  nozzle  exit  plane. 

Profile  measurements  of  the  mean  velocity  and  longitudinal  turbulence 
intensity  in  the  turbulent  boundary  layer  developing  at  20  m  s  indicate  normal, 
spanwise  uniform  behavior.  For  this  qualification  test,  and  all  results  which  follow,  the 


boundary  layer  was  tripped  near  the  exit  of  the  nozzle  with  a  1.5  mm  high  strip  of 
tape.  Total  pressure  measurements  along  the  test  section  surface  at  the  nozzle  exit 
were  uniform  within  0.5%  indicating  spanwise  uniform  skin  friction. 

Freestream  turbulence  intensity  was  measure  to  be  0.00085  (8.5  one  - 
hundredths  of  one  percent  or  .085  precent)  at  20  m  s  increasing  to  0,00095  at  30  m.s. 

B.  INJECTION  SYSTEM 

Ordinarily,  the  injection  system  provides  film  coolant  at  temperatures  above 
ambient.  The  freestream  air  is  at  ambient  temperature;  therefore,  the  heat  transfer 
would  then  be  in  a  direction  opposite  to  that  which  occurs  in  gas  turbines.  The 
coolant  is  injected  into  the  boundary  layer  through  a  single  row  of  injection  holes.  The 
injection  holes  are  scaled  relative  to  boundary  layer  thickness  to  be  similar  to  those 
used  in  current  turbine  blade  design. 

For  the  present  tests,  all  injected  air  was  at  ambient  temperatures. 

Injection  system  air  is  provided  by  an  71TD  Ingersoll-Rand  air  compressor.  The 
air  is  discharged  from  the  compressor  into  three  large  storage  tanks.  As  the  schematic. 
Figure  2.2,  shows  the  air  flows  from  the  storage  tank  through  an  adjustable  regulator, 
a  cut-off  valve,  moisture  separator,  flow  regulator,  a  Fisher  and  Portor  rotometer  (full 
scale  9.345E-3  m/s,  19.8  SCFM,  model  10A3565A).  The  rotometer  which  controls  the 
volumetric  flow  rate,  discharges  the  film  coolant  through  a  diffuser  and  into  the 
injection  heat  exchanger  and  plenum  chamber. 

The  heat  exchanger  and  plenum  chamber  shown  in  Figure  2.3  is  0.305  x  0.50S  x 
0.457  m  (23  x  20  x  18  in.)  and  is  constructed  of  1.27  cm  (1/2  in.)  plexiglass.  Injection 
air  flows  over  three  metal  plates  0.381  x  0.508  m  (15  x  20  in.).  The  two  lower  plates 
are  covered  by  silicon  rubber  heaters,  0.381  x  0.483  m  (15  x  19  in.),  rated  at  120  volts. 
The  heaters  are  controlled  through  a  type  136  Powerstat  variable  autotransformer. 
The  top  surface  of  the  chamber  contains  13  plexiglass  injection  tubes  each  being  8  cm 
(3.15  in.)  in  length  and  with  an  inner  diameter  of  0.95  cm  (3  8  in.).  This  corresponds 
to  a  length  to  diameter  ratio  of  8.42.  The  13  injection  holes  discharge  the  coolant  into 
the  boundary  layer  at  a  30  degree  angle.  There  is  a  three-diameter  spanwise  spacing 
between  the  center  of  each  hole. 

Further  discussion  of  the  qualification  and  performance  of  the  injection  system 
can  be  found  in  [Ref.  10:  p.23]. 


C.  FIVE-HOLE  PRESSURE  PROBE 


Multi-hole  pressure  probes  are  invaluable  in  the  investigation  and  measurement 
of  complex,  three-dimensional  (lows.  In  particular,  the  five-hole  probe  is  well  suited 
for  measurement  of  three  mean  velocity  components  in  low  speed  incompressible  flows. 

The  five-hole  pressure  probe  used  to  measure  pressure  in  this  study  is 
manufactured  by  United  Sensors  and  Control  Corp.  (drawing  number 
DA-125-24-F-22-CD).  The  probe  shown  in  Figure  2.4  is  0.6096  m  (2  ft.)  in  overall 
length  with  a  probe  diameter  of  0.318  cm  (0.125  in.).  It  is  constructed  of  corrosion- 
resistant,  non-magnetic  stainless  steel. 

The  five  pressure  holes  are  arranged  in  two  different  planes  which  intersect  at  the 
mutual  hole,  Pj.  The  probe  tip  is  prismatic  in  geometry,  as  shown  in  Figure  2.4  The 
centrally  located  Pj  hole  is  normal  to  the  freestream.  The  pitch  plane  consists  of  Pj, 
P^,  and  P^  pressure  holes.  While  pressure  holes  Pj,  P2.  P3  constitute  the  yaw  plane. 
The  distance  seperating  P2  and  P3  is  0. 1 7S  cm  (0.070  in.),  P^  and  P5  are  0.155  cm 
(0.061  in.)  apart.  The  central  hole  Pj  is  0.648  cm  (0.255  in.)  from  the  bottom  of  the 
probe  tip.  These  distances  were  measured  using  a  micrometer. 

For  calibration,  the  probe  was  positioned  in  a  manual  traversing  unit 
manufactured  by  United  Sensor  and  Control  Corp.,  Figure  2.5.  The  unit  was  modified 
by  the  user  to  include  a  compass  rose  with  a  radius  of  15.24  cm  (6  in.)  and  a  range  of 
yaw  angles  from  -40  to  +40  degrees  (0.25  degree  accuracy).  Additionally,  the  manual 
traversing  unit  positions  the  probe  vertically  in  the  mainstream  0  -  30.48  cm  (0  -  12  in.) 
range.  The  manual  traversing  unit  is  mounted  on  top  of  a  spanwise,  horizontal  sled. 
The  sled  is  62.23  cm  (24.5  in.)  wide  and  is  designed  to  set  in  place  on  top  of  the  wind 
tunnel's  side  walls.  The  spanwise  mounting  sled  is  designed  such  that  the  probe  can  be 
positioned  10.16  cm  (4  in.)  either  side  of  centerline  in  increments  of  0.635  cm  (0.25  in.). 
In  addition  to  spanwise  positioning,  the  spanwise  sled  could  be  rotated  through  a  range 
of  pitch  angles  from  -15  to  +15  degrees  with  an  accuracy  of  0.5  degrees. 

After  the  completion  of  the  probe  calibration,  an  automated  traversing 
mechanism,  Figure  2.6,  was  used  for  probe  positioning  while  measuring  the  pressures  in 
the  experimental  test  cases.  The  probe  is  fixed  into  the  automated  traversing 
mechanism  in  a  position  of  zero  yaw.  The  traversing  mechanism  has  two  degrees  of 
movement  which  allows  a  thorough  measurement  of  the  flow  field  to  be  conducted. 
Both  the  spanwise  and  vertical  traversing  blocks  are  mounted  on  a  20-thread  per  inch 
drive  screw  and  two  ground  steel,  case-hardened  steel  guide  support  shafts.  Each  drive 


shaft  is  directly  coupled  to  a  SLO-SYN  type  MO92-FD310  stepping  motor.  The 
motors  are  controlled  by  a  MITAS  Two-Axis  Motion  Controller,  Figure  2.7  The 
stepping  motors  and  the  controller  are  manufactured  by  Superior  Electric  Company. 
The  controller  directs  the  movement  the  probe  in  both  the  spanwise  amd  vertical 
directions  The  MITAS  controller  comes  equipped  with  2K  bytes  of  memory  and  an 
MC68000,  16-bit  microprocessor  which  allows  the  user  to  program  the  start,  stop, 
duration,  speed,  acceleration  and  deceleration  of  the  stepping  motors. 

D.  DATA  ACQUISITION  SYSTEM 

The  data  acquisition  system,  shown  in  Figure  2.8  rapidly  acquires  the  voltages 
associated  with  each  pressure,  converts  each  voltage  to  pressure. 

1.  Transducers  and  Demodulators 

The  probe  is  connected  through  reinforced  plastic  tubing  to  five  Celesco  model 
LCVR  differential  pressure  transducers.  These  transducers  have  a  designed  pressure 
range  of  0  -20  cm  (0  -  7.85  in.)  water  and  produce  a  15  to  45  mV  volt  output  signal. 
The  transducer  output  signal  is  converted  to  a  proportional  DC  signal  by  Celesco  CD 
10D  carrier  demodulators.  Each  demodulator  has  a  maximum  frequency  response  of 
-3dB  at  500  Hz  and  a  maximum  out  put  noise  of  lOmV,  peak  to  peak.  Each 
transducer  carrier  demodulator  combination  was  calibrated  against  a  Meridian  1.27  cm 
(0.5  in.)  horizontal  manometer  with  an  accuracy  of  0.002  cm  (.005  in.)  of  water,  to  give 
typical  calibration  of  approximately  1.0  volt  per  inch  of  water  differential  pressure. 

2.  Computers  and  Hardware 

A  Hewlett-Packard  85  microcomputer  was  used  to  acquire  and  process  data 
for  the  calibration  of  the  pressure  probe.  Configured  with  64K  bytes  of  memory  and  a 
single  magnetic  tape  cartridge  drive,  the  HP-S5  was  used  to  collect,  store,  display,  and 
print  the  majority  of  the  data  required  during  the  course  of  probe  calibration. 

For  the  measurement  of  the  (low  field,  a  Hewlett-Packard  Series  300,  Model 
9836S  computer  was  dedicated  to  the  data  acquisition  process.  The  HP  9S36S  is 
equipped  with  a  MC68000,  8  MHz  16  32-bit  processor.  Dual  5-1  4  inch  floppy  disk 
drives,  and  1M  bytes  of  memory  .  A  UP  7470  two  pen  plotter  was  used  for  the  graphic 
representation  of  data. 

Each  transducer  carrier  demodulator  combination  is  connected  directly  to  a 
HP  3498A  extender  which  is  controlled  by  a  HP-3497A  data  acquisition  control  unit. 
The  HP-3497A  which  provides  precision  measurement  and  process  monitoring,  is 
equipped  with  analog  multiplexing  and  a  digital  voltmeter  with  IpV  sensitivity. 
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Six  software  programs  were  developed  for  use  during  the  thesis.  PROCAL 
was  developed  for  use  with  the  HP-85  and  was  used  during  the  probe  calibration  phase. 
PRSACQ,  VEL,  VELC,  PLOT,  VECTOR  were  developed  for  use  with  the  HP  9836S. 
PRSACQ  is  used  to  measure  the  pressures  in  the  various  flow  fields.  VEL  and  VELC 
are  used  to  compute  the  velocity  components.  PLOT  and  VECTOR  are  used  for  the 
plotting  of  results.  A  thorough  discussion  of  the  programs  requires  an  understanding 
of  the  calibration  procedures  and  the  velocity  measurement  techniques  as  discussed  in 
Chapter  Three. 
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111.  CALIBRATION  AND  MEASUREMENT  PROCEDURES 


Prior  to  using  the  five-hole  pressure  probe  to  measure  flow  velocities,  it  is 
necessary  to  calibrate  the  probe  to  determine  the  dependence  of  yaw,  pitch,  static,  and 
total  coefficients  on  yaw  and  pitch  angles. 

A.  COORDINATE  SYSTEM 

A  right  hand  coordinate  system,  Figure  3.1  was  established  for  use  throughout 
the  course  of  study.  The  X-axis  is  parallel  to  the  streamwise  direction  and  is  positive  in 
the  downstream  direction.  The  Y-axis  is  in  the  vertical  plane  and  is  positive  from  the 
wind  tunnel's  bottom  wall.  The  Z-axis  is  in  the  vertical  plane  and  is  positive  from  the 
wind  tunnel's  bottom  wall.  The  Z-axis  is  in  the  spanwise  direction.  The  origin  of  the 
coordinate  system  is  located  on  the  centerline  line  of  the  bottom  wall. 

Yaw,  p,  is  defined  as  rotation  about  the  Y-axis  and  was  arbitrarily  defined  as 
positive  when  the  direction  of  the  flow  resulted  in  pressure  P-j  being  greater  than  P->. 
This  condition  corresponds  to  positive  component  of  velocity  in  Z  direction.  The  pitch 
angle,  a,  is  defined  as  rotation  about  the  Z-axis  and  is  defined  as  positive  when  P^  is 
greater  than  P^.This  condition  results  in  a  positive  component  of  velocity  in  the  Y 
direction. 

B.  EXPERIMENTAL  PROCEDURES 

The  five-hole  pressure  probe  was  calibrated  using  the  method  described  by 
Treaster  and  Yocum,  [Ref.  12]. 

To  perform  the  calibration,  the  probe  is  placed  in  the  manual  traversing  unit 
which  is  mounted  on  top  of  the  spanwise  horizontal  sled  as  described  in  Chapter  Two. 
The  probe  is  normal  to  the  freestream  when  P>  is  equal  to  P^  which  gives  P  equal  to 
zero. 

The  probe  was  manually  fixed  at  a  predetermined  yaw  angle  and  then  rotated 
through  the  pitch  plane.  At  each  point,  the  data  acquisition  system  records  the  five 
pressures.  Four  pressure  coefficients  are  then  calculated  by  the  PROCAL  program. 

Calibration  of  the  probe  was  conducted  over  a  range  of  yaw  angles  from  -20 
degrees  to  +  20  degrees  in  four-degree  increments.  The  pitch  angle  was  varied  from  -15 
degrees  to  +15  degrees  in  five-degree  increments.  This  provided  a  cone  of  angles 
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which  was  sufficient  for  the  flows  to  be  studied  in  this  thesis.  Calibration  was 
performed  in  the  NPS  Shear  Layer  Tunnel  at  a  freestream  velocity  of  approximately  21 
ms  (68.89  ft,  sec).  The  probe  was  located  1.6  m  (5.3ft)  from  the  boundary  layer  trip. 
This  equates  to  a  Reynolds  number  of  1.97x10^  based  on  downstream  distance.  Static 
pressure  was  obtained  from  a  static  pressure  tap  on  the  tunnel  side  wall  and  total 
pressure  from  a  Kiel  probe  inserted  into  the  flow  through  the  top  wall.  Static  and  total 
pressures  are  measured  once  per  set  of  calibration  points.  Atmospheric  pressure  is 
used  as  the  reference  pressure. 

C.  CALIBRATION  COEFFICIENTS 

To  use  the  probe  for  measurement  of  complex  flows,  it  is  necessary  to  determine 
the  flow  angles,  a  and  p,  and  the  local  static  and  total  pressures.  This  can  be  done  by 
determining  four  non-dimensional  pressure  coefficients  over  a  range  of  angles  in  both 
the  yaw  and  pitch  planes.  The  four  calibration  coefficients  are  defined  as: 


CPyaw  =  (P2_P3^Pl  “P)  (eqn  3.1 ) 

CPpitch  "  (P4-P5)  (Pj  “P)  (eqn  3.2) 

CPtotal  *  <Pl  ”  Ptotal^Pi  ~  ^  (eMn  3  3) 

^"Pstatic  "  ^P ~  Pstatic^Pl  ~  P*  <e9n  3 

where 

P  *  ( P,  +  P3  +  P,j  +■  Pj )  4  ( eqn  3  5 ) 


To  be  of  value  in  measuring  a  How  field,  the  calibration  coefficients  must  be  .1 
function  of  flow  angle  only,  independent  of  velocity  and  repeatable  figure  V2  shows 
that  the  calibration  is  both  independent  of  velocity  and  repeatable  in  the  \aw  plane 


As  previously  stated,  probe  calibration  was  conducted  by  fixing  the  yaw  angle 
and  varying  the  pitch  angle.  Theoretically,  a  calibration  procedure  of  fixing  the  pitch 
angle  and  rotating  through  the  yaw  angles  should  provide  identical  results.  However, 
in  this  study  when  the  fixed  pitch  vary  yaw  method  was  attempted,  the  results  were  not 
identical.  The  trend  of  Cppj^  vs.  pitch  angle  showed  significant  scatter,  and  Cpt0[Jj 
vs.  pitch  angle  was  not  constant  for  each  yaw  angle. 

D.  RESULTS 

Figure  3.3  shows  the  variation  of  Cpvaw  with  the  yaw  angles  for  various  pitch 
angles.  The  response  is  nearly  linear  for  all  yaw  angles.  The  results  show  that  values 
of  Cpvaw  for  various  pitch  angles  collapse  on  top  of  each  other  for  yaw  angles  of -12 
degrees  to  +•  8  degrees.  This  means  that  Cpvaw  is  independent  of  pitch  angle  in  this 
range.  From  -23  degrees  to  -15  degrees  and  from  +8  degrees  to  -*-15  degrees,  there 
are  slight  variations  in  the  results  indicating  that  Cpv  ,  has  a  slight  dependency  on  the 
pitch  angle.  The  fact  that  yaw  is  only  slightly  dependent  on  the  pitch  angle  allows  data 
to  be  more  easily  processed  in  determination  of  (low  velocity. 

The  variation  of  CPpjtch  vs.  pitch  angle.  Figures  3.4  and  3.5.  show  that  this 
coefiicient  is  dependent  on  both  yaw  and  pitch  angles.  The  trend  of  UPp,^  vs.  pitch 
angle  is  generally  linear  but  there  are  variations  from  that  linearity  for  each  saw  angle, 
and  unlike  the  yaw  plane,  these  variations  are  not  restricted  to  any  particular  region. 

Figures  3  4  and  3.5  show  that  the  range  of  values  for  Uppltch  is  small  over  the 
range  of  a  shown  compared  to  the  variation  of  Up  v  with  f$  [Ref  12:  pp  2"-2S| 
attnbtues  this  to  the  types  of  surfaces  on  which  the  holes  in  the  pitch  and  yaw  plane 
are  connected  Large  yaw  angles  result  in  one  hole  being  nearly  aligned  with  the  (low. 
This  hole  senses  a  pressure  nearly  equal  to  the  total  pressure  of  the  frecstream.  I  he 
other  hole  is  then  blocked  from  the  freestream  and.  consequently,  reads  a  pressure 
much  less  than  the  freestream  static  pressure  The  holes  in  the  pitch  plane  ha\e  a 
diirerent  response  to  variation  in  the  pitch  angle  When  pitched,  one  hole  reads  a 
pressure  which  is  near  the  total  pressure  of  the  freestream.  hut  the  second  hole  senses  a 
pressure  which  is  greater  than  the  freestream  static  pressure  1  hus,  Cpv  is  a  much 
larger  number  than  Upp,^-  The  small  range  of  ( 'p  ^  increases  the  scatter  and 
uncertainty  of  pitch  angle  measurements. 

The  plot  of  CpstJtlc  vs.  pitch  angle.  I  igure  3  (>.  shows  that  <  PvtJtK  has  wc'ak 
dependency  on  yaw  and  pitch  angles. 


Cptotal  vs-  P^ch  angle,  Figure  3.7,  indicates  that  for  any  given  yaw  angle, 

Cptotai  is  constant  throughout  the  range  of  pitch  angles.  Here,  Cp{Qtaj  also  shows  a 
weak  dependency  on  yaw  and  pitch  angles. 

E.  INTERPOLATION  AND  APPLICATION 

After  the  probe  has  been  calibrated  and  the  operating  characteristics  of  the  probe 
are  known,  it  is  possible  to  determine  the  pitch  angle,  yaw  angle,  local  static  and  total 
pressures  for  any  (low  field, 

The  probe  is  positioned  normal  to  the  freestream,  and  at  any  location  in  the  flow 
field,  the  five  pressures  can  be  measured.  These  five  pressures  are  now  used  to 
calculate  the  experimental  or  local  coefficients  of  yaw.  pitch,  static  pressure,  and  total 
pressure. 

A  fifth  order  polynomial  was  fitted  to  the  average  values  of  Cpvavv.  The 
resulting  polynomial  computed  using  a  FORTRAN  program  based  on  the  least-squares 
method  is: 

Papp«  -.158-7.36(Ml)  +  O.I35(VI1)2  +  0.3<M(M1)3  +  0.009(M1)4-0.03I(M1)5  (eqn  3.5) 

Here.  Papp  is  the  approximate  yaw  angle,  and  M  j  is  the  local  coefficient  of  yaw. 

Knowing  the  approximate  yaw  angle  and  the  local  CPpjtch  a  computerized 
interpolation  is  performed  to  determine  the  pitch  angle.  Since  the  value  of  the  pitch 
angle  is  dependent  on  the  local  CPpltch  and  on  the  yaw  angle,  it  is  necessary  to 
perform  a  double  interpolation.  Referring  to  Figure  3.8,  the  two  yaw  angles,  Pj  and 
Pj  which  bracket  the  approximate  yaw  angle  are  first  determined.  Next,  the  local 
Cppuch  M->a  and  is  located  between  known  values  of  Cppj^  from  the 

calibration  data  for  each  yaw  angle.  In  Figure  3.8,  these  values  are  designated  C|| 

Cj2  C->,  and  C22  which  correspond  to  pitch  angles  from  the  calibration  data  of  Oj  | 

Op  a-»j  and  The  bracketing  pitch  angle  a  j  f  is  determined  by  the  following 

interpolation 

(C)!  -  Vi2->  (CU  “  C12>“<ail  _air)  (an -ai2)  {e^n  3  6) 

which  rearranged  gives 

alf  “  a„  “(all  “al2K(Cu  -  Vi2a)  <<(  t,  -C|2>> 
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(eqn  3.7) 


The  other  bracketing  pitch  angle  c^p  is  found  in  a  similar  manner  which  gives 


a2f  “  a21  ■  (a21  “  a22^(c21  ”  M2b)  <C21  “  C22»  (eCln  3  8) 

The  pitch  angle  for  the  flow  dp  is  found  through  a  second  interpolation  which  results  in 
the  relationship 

®if  —  (®if  ~  ®2f^Pl  —  Papp^Pl  ~  (e<ln  J-9) 

Because  Cpvaw  has  a  slight  dependency  on  pitch,  it  is  necessary  to  refine  the 
approximate  yaw  angle  once  the  local  pitch  angle  is  known.  This  is  done  with  a 
computerized  interpolation  routine  which  is  very  similar  to  that  used  for  the  pitch 
angle.  As  shown  in  Figure  3.9,  the  pitch  angles  from  the  calibration  data,  ct|  and  a-, 
which  bracket  the  local  pitch  angle,  tip  are  first  determined.  The  local  Cpvau.  M.u  is 
then  located  between  the  values  of  Cpvavv  from  the  calibration  data  for  a  j  This  results 
in  C^j  and  Ctj  which  correspond  to  yaw  angles  pjj  and  Pp.  Then,  by  linear 
interpolation  Pjpis  found 

<C3l  “  Mib>  <C3i  ~C2|)~(((Pn  -Pif)  (Pi,  -  Pt2»  (eqn  3. 10) 

or 

Pi f *  Pi i  ~ lP[ [  —  P  12^^31  ~  ^31  —  ^"21^  (CflU  3.11) 

The  same  interpolation  is  done  for  a  2  which  gives 

P;f*p2l-<p2l"P22)«C32"Mla)(C32_C22))  ,c^n  U2) 

p->p  and  P2p  are  the  yaw  angles  which  bracket  the  yaw  angle  of  the  flow.  The 
flows  yaw  angle,  pp  can  be  found  by  a  second  interpolation  which  results  in  the 
relationship 

Pr-p2r-<P2r-pir)<<a,-af)  (a, -a,  >>  (eqn  3.13) 
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Since  local  total  and  static  pressures  are  not  measured  at  each  probe  location,  it 
is  necessary’  to  perform  another  double  interpolation  to  compute  their  values.  From 
the  calibration  data,  the  values  of  Cptotaj  for  each  of  the  bracketing  yaw  angles  and 
pitch  angles  are  known.  In  Figure  3.10,  these  points  are  designated  Cjj.  Cj2>  C-»j, 
and  C22-  An  interpolation  is  performed  to  determine  the  Cptola|  corresponding  to  the 
local  pitch  angle.  These  points  are  designated,  Cjf  and  C2P  The  local  Cptotaj,  M^,  is 
found  using  the  relationship 

(eqn  3.14) 

which  when  rearranged  as  shown  below  gives  the  local  Cptotaj  for  the  flow 

M3  =  Clf-  (C,f-  C2fX(Pj  "  Pf)  (Pj  -  P2))  (eqn  3. 15) 

The  local  Cpstatjc  is  found  in  a  similar  manner.  Referring  to  Figure  3.11,  the 
relationship  for  determining  the  local  Cp  statjc  is  then  given  by 

“  C3f  (C,f  C4fK(p!  ~Pf)  (Pj  -P2))  (eqn  3.16) 


F.  VELOCITY  COMPONENT  DETERMINATION 

Once  the  values  of  the  local  Cptotaj  and  Cpstatjc  have  been  determined,  the  local 
total  and  static  pressure  can  be  calculated  and,  subsequently,  the  total  velocity  at  the 
probe  tip  can  be  determined. 

The  defining  relationships  for  CptQtaj  and  Cpslatjc  can  be  rearranged  to 
determine  Ptolal  and  P$tatlc: 


^total=^i  (CPtotal^i  ^ 


(eqn  3.17) 


^static  **  ^Pstatic^^i  ^ 


(eqn  3.18) 


By  using  Bernoulli's  equation,  the  magnitude  of  the  local  total  velocity  is: 
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^  'v^*>total~"  ^static)  P 


(eqn  3.19) 


The  three  components  of  velocity  can  now  be  determined  using  the  total  velocity 
vector  and  the  local  pitch  and  yaw  angles.  Referring  to  Figure  3.1,  these  velocity 
components  are  given  by, 


V  =  VcosacosP 


V  =  Vsina 

> 


Vz  =  Vcosasinp 


G.  SPATIAL  RESOLUTION  CORRELATION 


(eqn  3.20) 


(eqn  3.21) 


(eqn  3.22) 


The  calculation  of  Vv  may  be  influenced  by  the  local  total  velocity  gradient. 
Corrections  for  this  effect  may  be  made  using  the  following  relationship 


Vyf=  vyo  +  (^U'5xK,y) 


(eqn  3.23) 


Here  V  is  the  uncorrected  value  ofV v  and  Vvf  is  the  value  of  Vv  corrected  for 
spatial  resolution.  The  value  of  ly  used  was  slightly  greater  than  the  distance  between 
P^j  and  (0.155  cm):  a  value  of  0.200  cm  gave  constant  Vv  through  the  two- 
dimensional  boundary  layer. 

H.  SOFTWARE 

Six  programs  were  developed  for  use  during  this  study.  They  are  PROCAL, 
PRSACQ,  VEL  VELC,  PLOT,  and  VECTOR.  Each  program  is  written  in  BASIC. 
PROCAL  was  used  with  the  HP-S5,  all  others  were  written  for  the  HP-9S36S.  All  of 
the  programs  are  listed  in  Appendix  C. 

PROCAL  is  a  BASIC  program  used  for  the  calibration  of  the  pressure  probe. 
The  program  begins  by  computing  the  correction  factor  for  random  noise  associated 
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with  each  transducer.  The  user  is  then  prompted  to  manually  calibrate  each  transducer 
against  a  horizontal  monometer.  Static  pressures  are  input  utilizing  a  static  pressure 
tap  on  the  side  wall  and  one  of  the  transducers.  Total  pressure  of  the  freestream  is 
input  using  a  Keil  probe  inserted  into  freestream  and  a  transducer.  The  user  is  next 
prompted  to  input  the  ambient  pressure  in  inches  of  mercury.  After  the  user  has 
positioned  the  probe  at  the  desired  angles  of  yaw  and  pitch,  those  angles  are  input  into 
the  program. 

The  computer  then  acquires  the  five  voltages  from  the  data  acquisition  system, 
converts  each  voltage  to  a  pressure  in  inches  of  water,  and  then  calculates  the  yaw. 
pitch,  total,  and  static  coefficients  of  pressure.  Finally,  the  yaw  angle,  pitch  angle,  and 
the  four  coefficients  are  stored  on  a  separate  file  and  printed  out  utilzing  the  HP-85's 
internal  printer. 

PRSACQ  was  used  to  acquire  the  pressures  during  the  experiment.  PRSACQ 
begins  by  prompting  the  user  for  the  number  spanwise  and  vertical  data  points  and  the 
resolution.  A  matrix  of  data  point  location  is  then  computed.  Next,  the  transducers 
are  corrected  for  noise  and  calibrated  against  the  manometer.  Freestream  static  and 
total  pressures  are  measured  and  ambient  conditions  are  input.  The  program  enters  a 
loop  which  samples  each  pressure  ten  times  per  probe  location.  The  local  Cpvaw  and 
Cpp^ch  are  computed.  Probe  position,  Cpyaw  and  CPpjtch,  Pj  and  Ptotai  are  stored 
in  a  matrix.  At  the  end  of  the  data  collection  run,  these  values  are  read  into  a  data  file 
on  a  floppy  disk. 

VEL  is  used  to  process  the  raw  data  acquired  by  PRSACQ.  The  data  is  first  read 
into  computer  memory.  The  program  computes  an  approximate  Cpvavv  using  a 
polynomial  fit.  Next,  double  interpolation  subroutine  is  used  to  comput  local  Cppj^. 
A  second  interpolation  subroutine  is  used  to  refine  the  value  of  Cpvaw.  The  values  of 
the  Cppjtc^  and  Cpyaw  are  used  in  a  third  interpolation  subroutine  to  compute  the 
local  Cpstatjc  and  Cptotaj.  The  total  pressure  is  found  from  the  definition  of  CptQtaj. 
The  velocity  at  the  probe  tip  is  computed,  and  the  x,  v,  and  z  components  of  velocity 
are  determined.  Probe  position  (y  and  z  coordinates),  total  pressure,  total  velocity,  and 
V  Vv  and  V  are  stored  in  a  matrix  and  then  read  into  a  data  file. 

VELC  corrects  the  Vv  component  for  spatial  resolution.  PLOT  is  used  to 
generate  graphs  of  streamwise  velocity,  total  velocity,  and  total  pressure.  VECTOR  is 
used  to  plot  the  secondary  flow  vectors. 


IV.  EXPERIMENTAL  RESULTS 


The  study  was  conducted  in  three  parts.  The  first  was  the  measurement  of 
baseline  data  consisting  of  boundary  layer  profiles  in  a  two-dimensional  mean  How 
field.  The  second  was  the  measurement  of  the  boundary  layer  characteristics  at  a 
freestream  velocity  of  15  m,  s  with  1)  film  cooling  only.  2)  with  embedded  vortex  only, 
and  3)  with  both  film  cooling  and  embedded  vortex.  The  third  part  was  the 
measurement  of  boundary  layer  with  vortex  at  20  m  s,  and  measurement  of  boundary 
layer  with  vortex  and  film  cooling  at  a  freestream  velocity  of  1 1  m/s. 

A.  BASELINE  RESULTS 

The  baseline  boundary  layer  profiles  were  conducted  at  a  freestream  velocity  of 
22  nvs.  For  these  tests,  the  top  wall  of  the  wind  tunnel  was  adjusted  at  20  m  s 
freestream  velocity  so  that  a  zero  pressure  gradient  existed  within  0.15  mm  water  along 
the  length  of  the  test  section.  Profiles  were  taken  at  three  spanwise  locations  z=  +  2.54 
cm.  z  =  0  cm,  and  z  =  -2.54  cm.  Figure  4.1  shows  measurements  of  the  streamwise 
velocity,  Vx<  Figure  4.2  shows  measurements  of  Vr  Figures  4.3  and  4.4  show 
mesurements  of  the  normal  velocity,  Vy.  The  results  of  the  baseline  measurements 
indicate  behavior  expected  of  a  2-D  turbulent  boundary  layer,  since  the  figures  show 
mean  profiles  to  be  spanwise  uniform  for  all  three  velocity  components.  Figures  4.2 
and  4.3  show  that  the  Vz  and  the  corrected  Vy  components  are  small  and  nearly  zero 
as  would  be  expected.  Figure  4.4  is  the  plot  of  Vy  uncorrected  for  spatial  resolution. 

B.  15M/S  FREESTREAM  VELOCITY  RESULTS 

The  investigation  of  the  flow  field  was  conducted  by  using  the  five-hole  pressure 
probe  to  measure  pressures  at  800  points  in  a  spanwise  plane.  Data  was  taken  at  20 
different  vertical  locations,  each  having  40  spanwise  locations.  All  15  m  s 
measurements  were  taken  at  a  location  of  1.49  meters  from  the  boundary  layer  trip  or 
0.39  meters  from  the  heat  transfer  plate  leading  edge.  The  film  cooling  cases  were 
conducted  with  injection  air  at  75%  of  full  scale  on  the  rotometer  which  corresponds 
to  a  blowing  ratio  (ratio  of  coolant  to  mass  lluxes)  of  0.50. 


1.  Boundary  layer  with  film  cooling 

The  results  of  boundary  layer  with  film  cooling  only  and  no  vortex  are  shown 
in  Figures  4.5  to  4.8.  Figure  4.5  shows  contours  of  the  streamwise  velocity,  Figure  4.6 
shows  total  velocity  contours.  Figure  4.7  is  the  plot  for  total  pressure,  and  Figure  4.8 
shows  the  secondary  flow  vectors.  Away  from  the  wall,  outside  the  boundary  layer,  the 
first  three  figures  show  spatially  uniform  behavior  for  Vx,  V,  and  PQ.  Secondary  flow 
vectors  are  very  small  everywhere  in  Figure  4.8.  Near  the  wall  deficits  of  Vx,  V,  and 
PQ  correspond  to  locations  of  the  film  cooling  jets  which  are  located  about  every  3.0 
cm  from  tunnel  centerline. 

2.  Boundary'  layer  with  vortex 

The  vortex  is  generated  by  using  a  half-delta  wing  which  is  3.0  cm  high  with 
7.5  cm  chordand  an  angle  of  18°.  It  is  identical  to  vortex  generator  #  2  described  by 
Joseph,  [Ref.  10:  p.76].  The  vortex  generator  was  located  at  a  spanwise  location  of 
z  =  4.79  cm  (note  that  the  direction  of  +z  in  [Ref.  10:  p.75j  is  reversed).  Figure  4.9 
show's  the  streamwise  velocity  results,  Figure  4.10  is  a  plot  of  the  total  velocity 
contours,  Figure  4.11  is  a  contour  plot  of  total  pressure,  and  Figure  4.12  shows  the 
secondary  flow  vectors.  The  contour  plots  for  VX,V,  and  PQ  show  significant  deficits 
caused  by  the  generator  wake  which  is  rolled  up  with  the  vortex.  The  center  of  the 
wake  is  located  at  y=»3.3  cm  and  z  =  -3.05  cm.  Figure  4.12  show’s  that  the  vortex 
center  is  located  near  the  same  location,  with  overall  characteristics  similar  to  a 
Rankine  vortex.  Figure  4.13  shows  the  streamwise  vorticity  contours,  where  the 
vorticity  is  calculated  using 

(Of  =  (5Vz/dy)-  (dvydz)  (eqn  4.1) 

The  vorticity  is  largest  near  the  vortex  center  as  expected.  The  circulation  of  the 
vortex  is  estimated  to  be  0.2708  m  /s  using  the  equation 

r  =  jo)xdA.  (eqn  4.2) 

3.  Boundary  layer  with  vortex  and  film  cooling 

Figure  4.14  shows  the  streamwise  velocity  contours  for  this  case.  Figure  4.15  is 
the  total  velocity  contour  plots,  Figure  4.16  shows  the  results  for  total  pressure. 
Figure  4.17  shows  the  secondary  flow  vectors,  and  Figure  4. IS  is  a  contour  plot  of 


vorticity.  Figure  4.15,  4-16,  and  4-17  show  that  the  deficits  for  Vx,  V,  and  PQ  from  the 
film  cooling  are  no  longer  present  near  the  vortex.  This  result  shows  that  the  effects  of 
film  cooling  are  decimated  by  the  vortex.  The  effect  of  the  vortex  is  particularly 
pronounced  near  its  downwash  side.  As  for  the  previous  case,  Figure  4.18  shows  that 
vorticity  is  again  highest  near  the  vortex  center.  Figure  4.19  from  [Ref.  10:  p.99j  shows 
that  the  results  of  this  case  are  consistent  with  previous  work  by  Joseph  which  shows 
high  heat  transfer  rates  on  the  downwash  side  of  the  vortex  and  low  heat  transfer  rates 
on  the  upwash  side  of  the  vortex.  The  circulation  for  this  case  is  estimated  to  be 
0.2708  m^/s. 

C.  20  M/S  AND  10  M/S  FREESTREAM  VELOCITY  RESULTS 

Figure  4.20  is  a  plot  of  streamwise  velocity  for  the  case  of  an  embedded  vortex 
with  no  film  cooling  with  a  freestream  velocity  of  20  m/s.  Figure  4.21  is  a  plot  of  the 
secondary  flow  vectors  for  the  same  case.  Figure  4.22  and  Figure  4.23  are  for  a 
boundary  layer  with  embedded  vortex  and  film  cooling  at  11  m/s  freestream  velocity. 
The  blowing  ratio  for  the  20  m/s  case  is  0.38,  and  the  blowing  ratio  at  11  ms  was  0.68. 
The  results  for  both  cases  show  trends  which  are  similar  to  those  discussed  above.  The 
results  at  a  freestream  velocity  of  11  m  s  are  believed  to  be  less  reliable  because  of 
disturbances  at  the  inlet  of  the  tunnel  which  propagated  to  the  test  section  during  the 
time  period  the  data  was  acquired. 


V.  SUMMARY  AND  CONCLUSIONS 


Measurements  are  presented  of  boundary  layers  with  embedded  vortices  and  with 
film  cooling  for  freestream  velocities  of  15,  and  11  m,'s.  Measurements  of  a  boundary 
layer  with  embedded  vortex  and  and  no  film  cooling,  and  of  a  boundary  layer  with  film 
cooling  but  no  vortex  are  presented  for  freestream  velocity  of  15  m  s.  Plots  of  total 
velocity,  V,  streamwise  velocity,  Vx,  secondary  flow  vectors,  total  pressure.  PQ,  and 
streamwise  vorticity  are  presented  for  many  of  these  test  conditions. 

The  results  show  that  the  embedded  vortices  completely  dominate  the  flow  field 
in  boundary  layers  with  film  cooling.  This  is  indicated  from  the  plots  of  V,  Vv  and  PQ 
which  show  the  effects  of  film  cooling  to  be  completely  decimated  in  the  vicinity  of  the 
vortex.  This  result  is  consistent  with  the  heat  tranfer  results  of  Joseph,  [Ref.  10:  p.54], 
which  shows  a  localized  hot  spot  at  the  wall  near  the  same  location.  Future  film 
cooling  injection-hole  arrangements  in  turbine  blades  must  be  designed  to  compensate 
for  hot  spots  due  to  the  vortices. 

In  order  to  conduct  this  study,  a  five-hole  pressure  probe  was  calibrated  for  pitch 
and  yaw.  The  probe  was  then  used  to  measure  five  pressures  associated  with  the  flow. 
From  these  pressures,  total  velocity  and  the  x,  y,  and  z  components  of  velocity  were 
determined. 

A  boundary  layer  profile  was  conducted  to  verify  the  calibration  of  the  pressure 
probe,  measurement  procedures,  and  velocity  computations.  The  results  show  expected 
boundary  layer  behavior  with  a  small  Vy  and  Vz  component. 

It  is  recommended  that  flow  visulization  study  of  the  interaction  between  the 
vortex  and  film  cooling  be  conducted  to  enhance  the  understanding  of  this  complex 
phenomena. 
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F  igure  1.2  Smoke  and  oil  visulization  lor  rotor  blade 
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ure  3.2  Repeatability  of  calibration  results  at  different 
frcestream  velocities  lor  Cp,.„,.r  vs.  vaw  ansles. 
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Figure  3.7  Probe  calibration,  CptoU|  vs.  pitch  angles,  frccstrcam  velocity  of  31  m  s. 
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Figure  4.7 


Total  pressure  for  boundary  laver  without  embedded  vortex, 
75%  film  cooling,  I reestream  velocity  15  m/s. 
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SECONDARY  FLOW  VECTORS 


Figure  4.12  Secondary  flow  vectors  lor  boundarv  layer  with  embedded  vortex, 
without  film  coolin'.’.  Ircestrcam  velocity  15  m  s. 
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Figure  4.13  Vorticitv  contours  for  boundarv  layer  with, embedded  vortex 
without  'film  cooling,  freestream  velocity  lo  m, s. 
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Figure  4.14  Streamwise  velocitv  for  boundarv  layer  with  embedded  vortex 
and  film  cooling,'  freestream  velocity  15  m/s. 
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Figure  4.16  Total  pressure  for  boundary  layer  with  embedded  vortex 
and  film  cooling,  freestream  velocity  15  m/s. 


TO 


WWW 


BEDDED  VORTEX 
’5*  FILM  COOL] 


TRERMWISE  VORTICITY 


Figure  4.18  Vorticity  contours  for  boundarv  laver  with  embedded  vortex 
and  film  cooling,  freestream  velocity  15  m/s. 
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Figure  4.19  Stanton  number  ratio  for  vortex  at  z=  -4.79cm  with  him  cooling 
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Figure  4.20  Streamwise  velocity  for  embedded  vortex 
without  film  cooling  ,  freestrcam  velocity  20  m/s. 
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Figure  4.23  Secondary  How  vectors  for  embedded  vortex 
with  film  cooling.  Ireestrcant  velocity  1 1  m  s. 


*  l 

» 


74 


'Va'a-Va'aVCv’a'a’a  Lf.-Cv  .  I'CvVfc 


-V-" 


APPENDIX  C 
SOFTWARE 


•  f?c\M  PROGRAM  PPQChl 
'•O  REM  THIS  PROGRAM  COMPUTES 
O.-i  R£M  rHE  COEFFICIENTS  OF  VAN 
40  REM  AND  PITCH  FOR  THE  FIVE 
50  REM  HOLE  PRESSURE  PROBE 
60  REM 

TO  REM  DAVID  EVANS  NOV  1986 
80  REM 

38  REM  VARIABLE  NAMES 

100  REM 

110  REM  E<!>  IS  THE  VOLTAGE 
ISO  REM  READ  FROM  THE  DATA 
138  REM  ACOU 1ST  ION  SYSTEM 
140  REM 

1  5 0i  REM  A ■*  I  >  IS  THE  C ON V E R S  I  0 N 
ICO  REM  FROM  VOLTAGE  TO  INCHES 
155  REM  OF  HATER 
170  REM 

180  REM  P < I >  WHERE  1  =  1  TO  5 
198  REM  ARE  ThE  PRESSURES 
288  REM  FROM  THE  5  HOLES 
218  REM  OF  THE  PRESSURE  PROBE 
228  REM 

230  REM  PG  =  STAT  I C  PRESSURE 
248  REM  P7=T0TAL  PRESSURE 
258  REM  P  ,6>=AVEPAGE  PRESSURE 
268  REM  D 1 =  P2-P3 
278  REM  D2=P4-P5 
230  REM  03=P1  -P<.6> 

298  REM  Cl=  YAW  COEFFICIENT 
300  REM  C  2  =  PITCH  COEFFICIENT 
310  REM  C3=C0EFF ICIENT  OF  TOTAL 
320  REM  PRESSURE 

338  REM  C4=C0EFF ICIENT  OF  STATIC 

PRESSURE 

348  REM 

350  DIM  Cl  <  100>  .•  C2<  100,' ,  C3UU0  > 
360  DIM  C4  < 1 00 >  >  X 1 < 1 00 )» Y(1 00 ) 
370  REM 

330  DISP  "ENTER  DATE , < MMDDYY > " 
390  INPUT  Mi 

400  DISP  "ENTER  T I  ME  ,  < HHMM  > " 

410  INPUT  M2 
420  REM 

430  PRINT  Ml , , M2 
440  REM 

450  CREATE  '* CALDA 1  "  ,  2 
460  ASSIGN.**  1  TO  "CALDA  l" 

470  '  PRINT#  1  •  Ml, M2 


REM  TH I S  SECTION  COMPUTES 
PEM  THE  CONVERSION  FACTOR 
PEM  FOR  VOLTAGE  TO  INCHES 
PEM  OF  WATER  THE  USER  MUST 
PEM  INPUT  THE  MANOMETER 
FEM 

REM  THIS  SECTION  COMPUTES 
615  REM  _THE  ZERO  OFFSET ‘COR 
616  PEM  RECTION  FOR  THE 
TRANS- 

PE  M  THE  ZERO  OFFSET  COR- 
622  REM  RECTION  FOR  THE  TRA 
NS-  623  REM 
REM 

DISP  -COMPUTE  ZERO  OFFSET  " 
DISP  -CORRECTION  *' 

CUSP  '■  1  !  1  DISCONNECT  ALL 

TUBING  TO  THE  TRA 
OUCERS  ‘  M  '• 


NS  OUCERS*  Mn 

DISP  *'  HIT  CONT  WHEN  READY" 
PAUSE 
CLEAR 

FOR  1=20*3  TO  284 
Z2  =  0 

FOR  J=1  TO  59 

OUTPUT  709  “A  I  “  ;  I  ,  "VTl  " 

ENTER  709  ;  X2 

Z2=Z2+X2 

NEXT  J 

Z3  =  Z2-''50 

Z4 » I  - 1 99  >  =23 

NEXT  I 

REM 

REM 

REM  *******LOOP#l ********* 
REM 

FOR  1=208  TO  204 
11=1-199 

0I3P  "TRANSDUCER  CALIBRATION 
ENTER  THE  MANOMETE 
R  PRESSURE  UN 


R  PRESSURE 

H  2  0  >  " 

DISP  "ENTER  THE  PRESSURE" 
DISF  " TRANSDDUCER  HP  ".II 
INPUT  H2 
Z0  =  0 
REM 

REM  *****L00P#2****t****» 
FOR  J=i  TO  50 
OUTPUT  709  "  A  I"  ,  I  ,  "  V  T  1  " 
ENTER  70?  . 

Z0=Z0+X 
NEXT  j 

pem 


330  E 1 =  Z  0  5 0 

3  40  A  v.  1 -199  >=H2/ <E 1-Z4< I - 1 99 } 

350  PR  I  NT  "  A  "  .•  I  A  <  I  —  1  3  3  ) 

360  DI3F  "IF  A < I  »  IS  UNSRT 
PRESS  8  " 

370  INPUT  N2 

380  IF  N 2  =  8  THEN  GOTO  380 
330  ne:;t  t 
1080  REN 

1010  REM  ******* **************** 
1028  PEN 

1830  REM  ENTER  THE  VALUES  OF 
1048  REM  STATIC  AND  TOTAL 
1650  REM  PRESSURE. 

1068  REM 

1070  CliP  "INPUT  STATIC  PRESSURE 

■CINCHES  OF  H20  >  " 


1  0  S  0 

0I3P  ‘USE 

TRANSDUCER  HP  i" 

1  830 

OISP  "AND 

CHANNEL  NR  200“ 

x  1  00 

REM 

1110 

CUSP  "HIT 

CONTINUE  WHEN  " 

1  120 

DISP  "P  STATIC  PROPERLY" 

1  1  38 

OISP  " HOOKED  UP  " 

1  140 

PAUSE 

1  150 

CLEAR 

1  160 

86  =  0 

1  1  78 

FOP  J= 1  TO 

50 

1  180 

REM 

1  1  30 

OUTPUT  783 

• "AI " . 2^9 • " V T 1 " 

1  208 

ENTER  703 

;  X6 

1210 

B6=B6+X6 

1228 

NEXT  J 

1230 

P6=B6/50 

1248 

REM 

1258 

DISP  "INPUT  THE  TOTAL  PRES- 

0  >  " 

SURE , < INCHES  H2 

1260 

DISP  "  AND 

CHANNEL  NR  200" 

1278 

DISP  " 

ii 

1238 

OISP  "  HIT 

CONTINUE  WHEN" 

1230 

DISP  "  P  TOTAL  PROPERLY" 

1  388 

DISP  "  HOOKED  UP  " 

1310 

PAUSE 

l  320 

CLEAR 

1330 

REM 

1  348 

B7  =  0 

1350 

FOP  J=1  TO 

50 

t  368 

OUTPUT  703 

•  "  A I " . 200 •  " V T 1 " 

1378 

ENTER  703 

X  7 

77 


1380  B7=B7+X7 
1390  NEXT  j 
1400  P7=B7 / 50 

1419  REM 

1420  REM 

14  30  REM  ENTER  ft  M B I E  N  T  CONDITION 

r 

1440  REM 

1450  CUSP  "ENTER  PftMB I M  HG  >  " 
14t'0  INPUT  ftl 
1470  REM 
1480  T  =0 

1490  FOR  J=1  TO  50 

1500  OUTPUT  709  ■  "ft  I  “  109  ,  "UTI  " 

1510  ENTER  709  .  T1 

1520  T  =  T  +  T ! 

1530  NEXT  J 
1  5  4  O  T  2  =  T  /  5  O 
1550  GOSUB  3010 
1560  REM 

1570  REM  DEL  P  COMPUTED  FROM 

1580  REM  PT07AL-PSTATIC 

1590  REM 

1600  C=P7-P6 

1610  REM 

1620  REM 

1670  REM  CONVERSION  TO  SI  UNITS 
1640  ft  l=ft 1*3385  82 


1650 
1660 
1670 
1  680 
1690 
1  700 
1710 
1720 
1730 
17«0 
1750 
1760 
1770 
1  730 
1790 
1300 
1810 

1320 
18  70 

l  340 

1350 


C  =C  *  243  7 
T  4  =  77  +  273  1 5 
R  1  =ft  1  2S7*T 4  ) 

U1  =  <2*C-'R1  >  '  5 

REM 

REM 

REM  ENTER  THE  BIG  LOOP 
REM  FOR  REftDING  PRESSURES 
REM  ftNO  COMPUTING  COEFF 
REM 

REM  * *****L00P#3*********** 
FOR  P=l  TO  100 
K  1  -  K 

REM  ENTER  THE  VALUE  OF  THE 
REM  PITCH  ANGLE  AND  YAW 
REM 

DISP  "  ENTEF  THE  VALUE  OF 

THE  PITCH  ANGLE" 

INFUT  :<!•*- 
F  EM 

DISP  "ENTER  THE  '..'hLUE  OF" 
DISP  "THE  TAW  ANGLE  " 


i  c 


I860  INPUT  YCIO 

1370  DISP  “  PITCH  ANGLE  ENTERED 

I  S  H  ,  X  1  <  K  > 

1330  DISP  "YAM  ANGLE  ENTERED  IS 

’•  .  Y  <  K 

1390  D I 3P  " 

1 900  0 1 3 P  “  DO  YOU  WANT  TO  CHANG 
E  PITCH  OR  YAW  IF 

30  PRE  S3  1" 

1310  INPUT  N4 

1320  IF  N4= 1  THEN  1300 

1930  REM  THIS  LOOP  ACQUIRES  EACH 

1940  REM  PRESSURE  50  TIMES  AND 

1950  REM  TAKES  THE  AVERAGE. 

1 93 y  R E  M 

1370  DISP  "HIT  'CONTINUE'  FOP" 
1930  OISP  "DATA  ACQUISITION  " 
1390  PAUSE 
2000  CLEAR 
2010  REM 

2O20  REM  *t****L0GP#3********* 
2030  FOR  1=200  TO  204 
2040  Z1=0 

2050  REM  **t***LG0P#4********* 

2060  FOP  J  =  1  TO  50 

2070  OUTPUT  709  ; " A  I “ ; I , " V T 1 " 

2030  ENTER  70?  ,  X 

2O90  Z1=Z1+X 

2100  NEXT  J 


2300  REM 

2310  PRINT#  1  X 1  <  K  -•  V  <  K  >  ,  C  1  •.  K  > 

,  C2<K>  ,  C3<K>  ,  C4*.'»  '* 

2320  REM 
2338  REM 

2348  PRINT  "PITCH  ANGLE  IS%XhK 

t 

2350  PRINT  "YAW  ANGLE  IS",Y<K,- 
2360  PRINT 

2370  PRINT  “P1=“,P<1> 

23S0  PRINT  "P2= “ ■ P< 2 > 

2390  PRINT  ,,P3=“  P(3;- 
2400  PRINT  "P4=" ,P(4> 

2410  PRINT  "  P  5  =  "  ,  P  <  5  > 

2420  PRINT 
2430  PRINT 

2440  0  ISP  "LAST  POINTS?" 

2450  DISP  "IF  SO  PRESS  1" 

2460  INPUT  HI 

2470  IF  N 1 = 1  THEN  GOTO  2500 
2480  NE'iT  \ 

2490  REM  *********************** 
1705 

2500  REM 

2510  ASSIGN#  1  TO  * 

2520  REM 
2538  PRINT 

2540  PRINT  " **t*****r***** ****** 
t  ******* " 

2550  PRINT 
2560  PRINT 

2570  PRINT  "FIVE  HOLE  PRESSURE 

PROBE  CALI  BRAT 

ION  " 

2580  PRINT 
2590  PRINT 

2609  PRINT  "DATE  OF  RUN  IS ".Ml 

2610  PRINT  "TIME  OF  RUN  IS", M2 
2620  PRINT 

2630  PRINT 

2640  PRINT  "DENS  I  TV  <  KG y M3 ) " 

2650  PRINT  R 1 

2660  PRINT  MUELOCITY<M,'S>" 

2670  PRINT  U1 

2630  PRINT  "PPMBvN/M2 >  " 

2690  PRINT  A 1 

2700  PRINT  " T  RMB  v  C  > " 

2710  PRINT  T3 

2720  PRINT  "P  STATIC  <  I N  H20  ;.  " 
2730  PRINT  F'6 

2740  PRINT  "P  TOTAL  =  ‘  IN  H20  •  " 
2750  PRINT  P7 
2760  PRINT 

27  70  PRINT  "PIT  TAW  f p I  f 


2  7  3  9  R  £  M  X  X  X  X  X  t  L  0  0  P  #  5  X  X  X  X  X  X  X  X  X  X 
2790  PGR  1=1  TO  k 1 
2800  PRINT  US  I NG  2810  ;  X 1 < I > , V < 
I  ">  ,  C  i  <  I  >  >  C  2  (  I  > 

2310  I M ft G E  MOD..  2X,  MDO..  2X.MD.DDDD 
.2  X  .  f-lD  .DDDD 

2820  NEXT  I 

2870  REM  XXXXXXXXXXXXXXXXXXXXXX 
2  8 4  0  P PINT 
2350  PRINT 

2360  PRINT  "PIT  YfiW  CPS  CP 

•  II 

H 

2370  REM  ******L00P#6********.#* 
2830  FOR  1=1  TO  K 1 
2380  PRINT  USING  2900  ;  X  Id),  Y  < 
I  >  .•  C 3 <  I  >  .•  C 4<  I  ) 

2900  IMAGE  MOD, 2X, MOD, 2X, MD . ODOD 
,2  X , MD . DDOD 

2910  NEXT  I 

2920  REM  XXXXXXXXXXXXXXXXXXXXXXX 

X 

2970  PRINT 

2940  OISP  "  THANKS  FOR  USING  ME" 
2950  CUSP  “  HOPE  YOU  FOUND  WHAT" 
2980  DISP  "YOU  WERE  LOOKING  FOR" 
2970  DISP  "  REGARDS" 

2930  DISP  "  HP-35B" 

2990  END 
3000  REM 

7010  REM  VOLTAGE  TO  TEMPERATURE 
3020  REM  CONVERSION 
3030  REN 

3040  REM  V  < MV >  TO  T ( C> 

3050  REM 

306  0  E 1  =  T2  P 1 000 

3070  T3  =  26 . 573*E1-1  936879  PE 1  PE  1 
+  .  99785*E1  *E1*E1-  .  261277TE1 
PE  1  PE  1 PE 1 
3030  RETURN 


RfcM  r  RUuRfiM  t'ril’n'./d 

I'!  h  M  THIS  P  R  G  G  R  h  M  n  C  0  ‘J I R  E  S  THE  P  P  £  S  U  R  F  c-  F  R u " 
REM  ri  I-  IvL  HOi-L  FREE  jURE  Ph'jiE  ni.'U  •_'-0 *J T  £_ 
REM  OBTAINtU  Ph’tjn  Tiit.  F Rub RAM  "FR"J'_'AL  " 


iifJU.HR  / ,  FESRVAR  < ,  i  j 


50  REM 

OB T  A I hL u 

EC  REM 

70  REM 

AUTHOR  g 

CO  REM 

‘JO 

vMk  i 

1  I'iCj1; 

lie  \  L  tr': 

C  y  ■  )  * 

ijy  Kt.r  I 

ri '  I  )  !  b 

♦  T-j  pprt 

<  T  '  t 

14C  REM 

to  s 

i S3  REM 

PRESSURE 

*  -■  f  » !..  !  ! 

170  REM 

DS  “ST AT  I 

rj-J  he  un, 
MOLE  fPESZ-U* 

1  Al-'O  1  :  T  r  t 


REM  f'7  =  TOTAL.  PRESSURE  OP  THE  FPEESTRErVi 
i’LM  R b  /“AvEPASfc  HR L 5 EURE 
REM  D1 -F’E-PZ 


REM  D.3=P1  •  F(  u  ) 


Ci “CALIBRATED  YAW  CQEEFICI ENT 
C2=CAL IBRATEO  FITCH  COEFFICIENT 
C3=CAL ISRATEC  TOTAL  COEFFICIENT 
04=CAL iBRATEC  STATIC  COEF c I C IENT 


0  G  r! 

Mi “MEASURED  YAW  CO 

EFFICIl 

NT 

p h* 

UNMEASURED  PITCH 

C  u£FF  I C 

IE‘iI 

ppM 

M3-MEASURE0  TOTAL 

OOEFFiC 

I  ENT 

REM 

^“MEASURED  STATIC 

cceff  : 

LIEN 

r-'EM 

REM 

V -THE  L. C A L  Tu  r  n L 

uSLOcr 

V  \  '  C : 

REM 

VI “THE  :<  COMPONENT 

OF  vE'l. 

OCIT 

REM 

U2“THE  V  COMPONENT 

OF  '.’EL 

f  T 

REM 

V  Z =T HE  1  COMPONENT 

Or  VEL 

rt~  y  t 

REM 

REM 

PITCH  ANGLE 

FROM  T 

He  '.7 

REM 

Y“THE  yaw  angle  fr 

OM  THE 

C  ri  L  I 

REM 

PEM 

,-4  =  THb  PITCH  HNGLE 

OF  THE 

FLO' 

REM 

REM  VI “THE  VERTICAL  LOLA T ICN  jr  TnE  3R2?£ 
REM  Z0=THE  SPANWISE  L  Of.fi  T I  Off  OF  THE  PRICE  T 
REM  Y2*THE  YAW  ANGLE  OF  THE  FlOW 
REM 

REM  A 1 'AMBIENT  PRESSURE 
6 EM  ^ABSOLUTE  AMBIENT  TEMPI  i-ELVlN  > 

::EM  F  1  “DENSITY  (KO-TV  3  ' 

,.-M  P-'T  L  C'C'TT'  Pc‘  THE  e?.,Ec"';  tc-rftH 


1580 

REM  EN 

0  UP  LOOP  FOR  CAL  I  £" A  1  I N>j  T  RmHSUU'.  FFS 

1510 

REM 

1528 

rem  enter  the  values  of  static  ano  total  pressure  of 

i  530 

PRINT 

"INPUT  STATIC  PRES'UKE  Clp  FREE  STREAM  ■  IN.HEG  ! " 

1540 

print 

"USE  TRANSDUCER  f\h  .  1  rtilD  CHANNEL  .  232 " 

1550 

PRINT 

1580 

PRINT 

"HIT  CONTINUE  UHtN  P  STATIC  ! -"  Fr-.:-Ep  L SE"  UP 

1570 

pause 

1 5c  0 

P0ST0 

1  5  70 

•  -  :  •  [■*  J  = 

i  rc  sc 

1500 

OUTPUT 

739 ; “ A I “ ; 203; "VT 1 " 

1610 

ENTER 

709 ; KB 

1S2S 

BS =86  + 

k£ 

IF  30 

NEXT  J 

1540 

REM  £N 

2*  L V F  F OR  fVJUUlF  INU  c-T,-P  IC  PAL;"  .  *JF  1  ;>5E 2 T pEa: 

1  c,  c  g 

P9=56/ 

53 

1668 

PRINT 

IS  70 

PRINT 

1530 

REM 

1690 

PRINT 

INPUT  TOTAL  F'RES cLfL  uF  ML  l_SCEc  i  FL**N . : 

1700 

PRINT 

"USE  TRANSDUCER  N«J .  1  AND  CHANNEL  Ka .  2?3" 

’ HE  pc-EESTPE«- 


1710  PRINT  "HIT  CONTINUE  WHEN  P  TOTAL  PROPERLY  55: VP' 

1720  PAUSE 

1730  57=8 

1740  FOR  J=i  TO  50 

1 753  OUTPUT  70S ;  “ A I  "  ; 200 :  " V I !  ' 

17S0  ENTER  709; X? 

>770  g?=e?t-:<7 

1730  M£>:T  J 

1730  REM  END  OF  LOOP  FOR  ACQUIRING  '"OTA!  r  P?  •-(-  ~-F  '  T  PR  Ac* 

1830  P7=B7.'S0 

1312  REN 

i 823  REM  ENTER  AMBIENT  CONDITIONS 
1333  INPUT  "ENTER  FAMB  « IN. OF  Ho  r  ,A1 
1840  REM 

1850  REM  LOOP  FOR  ACQUIRING  7PM”  OF  RREESTkE-M  "IP  IHERMOCOUF LE 
1360  T- 0 

1870  FOR  J=1  TO  50 

1880  OUTPUT  703;  "AI';103; -‘JTl  " 

1330  ENTER  703; Ti 
1900  T=T+T 1 
1910  NEXT  J 
1920  REM  END  OF  LOOP 
1930  T2=T/S0 

i940  REM  CONVERT  VOLTAGE  TG  'EMPEF A  TORE 
1950  E1=T2*1000 

1998  T  j ~2 B  .  b7 3 * E 1  - 1 . 3  jo 6  1 9 * E  2  >  .  N  2 .  - 5  *  1 1  '  .j  -  . 5  i  2  ~  * E 1  - 
1970  REM 

1330  REM  CONVERSION  TO  SI  UNIT. 

19  30  T  4  - T  3  +  2  7  3 . 15 
2330  A  1  ■ A ; *3385.52 


V*JV'j-V\>  V  _'J\  ‘j 


*.->*•  j.  ■>>: vCvS 


010 

R  i  = 

A I / .287 

* T4  ! 

320 

c = •: ; 

P7-P9 J* 

2  43.7 

030 

REM 

FREEST 

REAM  VElOC 

I  Tv 

040 

U3= 

>.  2  *  C  /  f?  1 

>  ■  _  c 

053 

REM 

060 

PEM 

BEFORE 

EACH  DATA 

Run 

ENSURE  TiiriV  '■ 

370 

print  “Balance  pc  an 

D  P3 

. 

350 

nn  r  i 

r  h  i  i 

NT  "HIT 

CONTINUE 

WHEN 

ALL  TUBING  Hr 

033 

P  r*«J 

SE 

100 

REM 

LOOP  F 

JR  SALANUT 

M3  Pi 

ROBE 

110 

RDF: 

1-231 

TG  232 

rcc/v 


i:a  Gi -o 

13?  FOR  j=i  TG  53 


\  •$.  0*1 

OUTPUT  70 

153 

ENTER  7S3 

160 

G l *5 i r :< 

170 

NE  <  T  j 

lb'0 

X— 31/50 

1 3  ‘-o 

E  ■.  I-i39;  = 

233 

P 1  I  -  i  33  ■  = 

2  1 0 

NEXT  I 

»■»  ;7» 

PRINT 

T‘  “?P 

PRINT 

240 

r  P I N  T  '■  p  0 

253 

pFlMT  " P5 

263 

INPUT  “DO 

270 

IF 

2  C><J 

!  <EM 

253 

REM  ENTER 

300 

PEM  AND  C 

331 

PF INTER  [ 

PPIfIT  y 

303 

PRINTER  I 

313 

REM 

”  •'/>. 
e..  Vj 

K  3  =  O 

321 

i.  2=M3*N3 

330 

FOR  K=1  7 

331 

!<g*K 

340 

REM 

390 

WAIT  10 

400 

REM  LOOP 

410 

FOR  1-230 

420 

01*0 

430 

REM 

440 

FOR  J-l  T 

450 

OUTPUT  70 

460 

ENTER  739 

473 

Gi =&i  +  X 

4  30 

NE  T  i 

490 

REM 

D  Uj 

•*  “ 1  ?  1  l  KJ 

?£=••  ,  P (  2  ) 


;?.e  COM^utt 

nt'JD  7140  FF 


m 


M2 


Pi 


!ENTbR  THE  DAS  AND  SAMPLt  EACH  PRESSURE  10  TIME 


I  .y  ;i. 


THE  VALUE 


86 


-  i  r  'i 


:si3 

Ei  1-199  ■  =  X 

e=:g 

F  •:  I  - 1 39  I  -  1  33  '  *6  :  I  - 

•  1  39  .'-04 

252 1 

IF  I =  2 00  THEN 

i  r*  •*. 

P!  ■  K  >--?■  1-193  ■ 

*-.  3  _  — 

IF 

_  D  jc1 

NcxT  1 

f.  '5  £  yj 

RE:’  END  Dr  LOOP  FOP  >vj 

DU i p I  NO 

-  Sb:j 

C  ’1 

.2  S  6  * 

h  EM  D-.-1P'..TE  The  AvERAG 

E  OF  P  - 

-> * — ■-  -» 

ncr> 

-  -  ■ 

2530 

■  'JP  i.-..  " 2  5 

2530 

'  J  b  -  P  '*  i  .‘-25 

■*•••/>  ••% 
-  S3  v  U 

me* r  i 

2  £  t  :7| 

Po :  *■'  '=Db  4 

:s:o 

-£71  Erie-  LOOP 

25.30 

c  Fri 

2S40 

HEM  COMPETE  TP£  COEFFt* 

Tic  r»i  t  s  r 

2E50 

FcH 

25S 1 

vl»Pi 2  )-P:  2 

•  obii 

D2»F\  4  )-P(5  > 

2-00 

D  3--p  (  1  )  — r  5 1  K  i 

2710 

MKK  -D1--03  ■  DOC* 

vr*;.y  c 

“>  ^  1  •  j\ 

M2'K  i-02/03  1  LOOP! 

.  PITCH 

“•  “7  T  > 

-  . 

printer  is  ?ai 

2  7  40 

PRINT  USING  2750  :VK 

:  y  • ,  r-i  i . 

27  =  3 

IMAGE  MOD .  CD  .  2:< ,  MCD .  DO . 

.jajidd. 

2763 
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APPENDIX  C 

UNCERTAINTY  ANALYSIS 


Uncertainty  analysis  was  performed  using  the  method  originally  attributed  to 
Kline  and  McCormick,.  (Ref.  13J.  Let  8f  be  the  uncertainty  in  the  result  and  6j  §2  • 
.  .,  5n  be  the  uncertainties  associated  with  each  independent  variable.  The  uncertainty 
in  the  result  can  be  expressed  as 

8r  -  H((aR^x.)x(8.))jl/2  (eqnC.l) 

To  determine  the  uncertainty  of  the  pitch  angle,  a,  a  straight  line  approximation 
was  made  for  a 

a  -  m  +  b(Cppitch)  (eqn  C.2) 

The  following  independent  variables  were  determined:  5m  -  ±1.97  •,  *  ± 

2.32«,  and  5^  -  ±  0.04637.  The  uncertainty  of  a,  was  determined  to  be  ±  2.36«. 
The  high  uncertainy  of  the  pitch  angle  is  caused  by  the  probe  being  highly  sensitive  in 
the  pitch  plane. 

The  uncertainty  of  the  yaw  angle.p,  is  determined  in  a  similar  manner.  Once 
again,  a  straight  line  approximation  is  made. 

P  *  m  +  b(Cpyaw)  (eqn  C.3) 

In  this  case,  the  uncertainties  of  the  independent  variables  were  found  to  be  6m 

-  ±  0.954,  5u  -  0.5924,  and  8r_  -  ±  0.058.  From  these  values  the 

o  '-Pvaw 

uncertainty  of  P  was  calculated  to  be  ±  f  .29  • 
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